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We start from an effective Hamiltonian for Ru ions in a square lattice, which includes the on- 
site interactions between t^g orbitals derived from Coulomb repulsion, and a tetragonal crystal- field 
splitting. Using perturbation theory in the hopping terms, we derive effective Hamiltonians to 
describe the Ru02 planes of RuSr2(Eu,Gd)Cu20g. For undoped planes (formal valence Ru +5 ), 
depending on the parameters we find three possible orderings of spin and orbitals, and construct a 
phase diagram. This allows to put constraints on the parameters based on experimental data. When 
electron doping consistent with the hole doping of the superconducting Ru02 planes is included, 
we obtain (for reasonable parameters) a double-exchange model with infinite antiferromagnetic 
coupling between itinerant electrons and localized spins. This model is equivalent to one used before 
[H. Aliaga and A. A. Aligia, Physica B 320, 34 (2002)] which consistently explains the seemingly 
contradictory magnetic properties of RuSr2(Eu,Gd)Cu208. 

PACS numbers: 75.10.-b, 75.10.Jm, 74.25. Ha, 74.72.Jt 



I. INTRODUCTION 

In recent years, there has been much interest in ruthcn- 
ates because of their interesting magnetic and super- 
conducting properties. For example, SrRuO*3 is a fer- 
romagnet that orders at Tm — 165 Kji Si'2Ru04 is an 
exotic p-wave superconductor with transition tempera- 
ture Ts = 1.5 Kj£ and Sr3Ru207 presents a metamag- 
netic transition— and non-Fermi liquid behavior^ A close 
relationship between ferromagnetic (FM) exchange and 
triplet p-wave superconductivity is expected in analogy 
with 3 He (Ref. |5|) or from bosonization studies in one 
dimension&i 

RuSr2(Eu,Gd)Cu208 has a magnetic transition at 
Tm ~ 133 K, and a superconducting transition at T<? ~ 
33 K for Eu or T s ~ 15 - 40 K for Gd (depending on the 
conditions of preparation and annealing) A 9 Below Tg, 
superconductivity coexists with magnetic order, which 
was first believed to be FM^iliLlLH since the magnetiza- 
tion shows a rapid increase with magnetic field for fields 
below 5 T, and the inverse magnetic susceptibility at high 
temperatures yields a positive Curie constant = 100±3 
However, neutron diffraction experiments found su- 
perlattice reflections consistent with an usual antiferro- 
magnetic (AF) order with nearest-neighbor spins antipar- 
allel in all three directions^ This seems difficult to rec- 
oncile with the above mentioned magnetic properties, in 
particular with a positive Curie constant. Nevertheless, 
a double-exchange model could consistently account for 
these observations 

The crystal structure of RuSr2(Eu,Gd)Cu20s is sim- 
ilar to that of YBa2Cu307, except that Y is re- 
placed by Eu or Gd, and the CuO chain layer is re- 
placed by a square planar Ru02 layer, with resulting 
tetragonal symmetry, except for small distortions typ- 
ical of perovskites^ The sequence of layers perpendic- 
ular to the tetragonal axis is Ru02/SrO/Cu02/(Eu or 
Gd)/Cu02/SrO. Several experiments, like muon spin 



rotation, 9 magnetization^*!^ electron paramagnetic reso- 
nance, and ferromagnetic resonance^ demonstrate that 
the development of superconductivity does nor affect the 
magnetic order. This suggests that, at least as a first ap- 
proximation, the superconducting Cu02 planes and the 
magnetic Ru02 planes behave as separate entities related 
only by charge transfer, as it happens with CuO planes 
and chains in YBa2Cu306+a;>^ Band structure calcula- 
tions are consistent with this pictured From what is 
known for YBa2Cu306+ic! 15 ' 16 a superconducting critical 
temperature Ts ~ 30-40 K suggests a doping of slightly 
less than 0.1 holes per Cu02 plane. This implies a dop- 
ing of < 0.2 electrons to the Ru02 planes with respect to 
the formal oxidation states Ru +5 and O -2 . Taking into 
account a certain degree of Ru-0 covalence, this doping 
is consistent with x-ray absorption near-edge structure 
(XANES) experiments which suggest a Ru valence near 
4.6^£ This situation is at variance with the compounds of 
the Ruddlesden-Popper series (Cu,Sr)„ + iRu„03 n +i (like 
those mentioned above) for which the formal oxidation 
state of Ru is Ru +4 . 

The main features of the puzzling magnetic behavior of 
Ru02 planes in RuSr2(Eu,Gd)Cu20s were explained in 
terms of a double-exchange model in which Ru +5 spins 
have a strong Hund coupling with a band of itinerant 
electrons^ Within this picture, the undoped system 
presents usual AF ordering. Additional electrons form 
FM polarons that tend to align easily in the direction 
of an applied magnetic field. Consequently, in spite of 
the AF order, the magnetic susceptibility at tempera- 
tures T > Tm can be well described by \ — C/{T - 6), 
with > 0, in agreement with experiment. 12 Further 
support to the double-exchange model is brought by 
the negative magnetoresistance above Tm, or below Tm 
for high enough magnetic fieldiiSiiS While this model 
has been successful in explaining several properties of 
manganitesj22i2i where the itinerant electrons are in 3d 
e„ orbitals, there is so far no justification for its appli- 
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cation to ruthenates, where the relevant orbitals are the 
4d t 2g oneS | n i 22 i 23 in which case crystal-field effects are 
expected to be more important, and correlations should 
be smaller due to the larger extent of the Ru 4d orbitals 
in comparison with Mn 3d ones. 

In this work, we start from an effective Hamiltonian 
H for the Ru 4d t 2g orbitals in a square lattice, after 
integrating out the O orbitals. H includes all atomic 
Coulomb interactions, and a tetragonal crystal-field split- 
ting A. We treat H in perturbation theory in the ef- 
fective Ru-Ru hopping. Since the effective parameters 
are difficult to estimate, the quantitative validity of this 
strong-coupling approach is difficult to address. How- 
ever, for reasonable parameters our results are consis- 
tent with experiment, and confirm results previously ob- 
tained with the double-exchange model. In the case of 
Sr 2 Ru0 4 , which has been studied in more detail, there 
is a wide range of proposed parameters, but it is clear 
that the correlations are significant, and the system is 
believed to be in the intermediate-coupling regime—^ 
Notice that approaches that neglect quantum fluctua- 
tions should assume smaller interactions to avoid mag- 
netic ordering in Sr2Ru04, while instead Ru02 planes in 
RuSr 2 (Eu,Gd)Cu 2 0s do order magnetically at Tm- 

The paper is organized as follows. In Section II we 
present the model and discuss its parameters. In Section 
III we describe the eigenstates an energies of the local 
Hamiltonian. In Section IV we derive effective Hamilto- 
nians that describe spin and orbital degrees of freedom in 
the undoped case, after integrating out the charge fluc- 
tuations. Section V contains the phase diagram for this 
case. In Section VI we discuss the effective Hamiltonians 



for the doped case, and their relation to the double ex- 
change model. Our results are summarized and discussed 
in Section VII. 



II. THE MODEL 

We start with an effective model for the 4d i 2g or- 
bitals of Ru ions in a square lattice. It can be derived 
from an appropriate multiband model for Ru0 2 planes 
by a canonical transformation eliminating Ru-0 hopping 
terms2& or by the cell perturbation method if Ru-0 co- 
valence were important The Hamiltonian is 



(f) 



where if) contains the local interaction terms at site i, 
Hq F is a tetragonal crystal-field splitting, and Hh con- 
tains the hopping terms which we restrict to nearest 
neighbors. Since if) contains only intra-site interactions, 
we assume for it the same form as for an isolated Ru ion, 
neglecting spin-orbit coupling. This form can be calcu- 
lated in a straightforward way using known methods of 
atomic physics. 29,30 Expanding the Coulomb interaction 
term e 2 /|i"i — r 2 | in spherical harmonics, all Coulomb in- 
tegrals can be expressed in terms of Slater parameters 
F , F-2, and F 4 (as done earlier— for e g orbitals). Here 
we write if) using the Kanamori parameters (which seem 
to be more popular in condensed mattes 3 ^ 3 ") for ti g or- 
bitals: U = F + 4F 2 + 36^4, J = J' = 3F 2 + 20F 4 , and 
U' = U- 2J. Then, 



J 



H} = U n ia ^n ia i + i ^2 {U'n iaa n ll3ljl + Jd\ aa d\ &a ,di aal d l j ia ) + J' 



4aT d L|^/31^/3T > 



where rii a(7 = d\ aa di acJl and d iaa creates an electron in the i 2g orbital a (xy, yz or zx) with spin a at site i. 

I 



(2) 



Choosing z as the tetragonal axis, we write the crystal 
field term in the form 



1) 



(3) 



in such a way that it changes sign under an electron-hole 
transformation. 

Denoting by 5 = ±x, ±y the four vectors that connect 
a site with its four nearest neighbors, the hopping term 
has the form 



-t^lidi+^zxA™* + 4+y,yz* d iyz<r + H.C.) 

- t 'J2 d l+S: Xy Axy° ■ ( 4 ) 



Ma- 



Notice that, since we neglect the distortions, electrons 
occupying zx {yz) orbitals do not hop in the y (x) di- 
rection due to the symmetry of the intermediate O 2p 
orbitalsi^&S 

While the parameters of H are difficult to estimate, we 
expect that the order of magnitude of t and t' is near I /4 
e y 14 i 24 i 25 i 26 i 28 (see also Section VII). Since the exchange 
interactions are not expected to be strongly screened in 
the solid, one may estimate J from atomic spectra— 
From the low lying levels of Ru + (with three holes in the 
4d shell), we obtain F 2 ~ 863 cm -1 and F4 ~ 78 cm -1 , 
leading to J ~ 0.5 eV. Optical experiments in Sr 2 Ru0 4 
suggest that U ~ 1.5 eV 3 - Notice that the expectation 
value of the Coulomb repulsion in any state with two 
electrons and total spin 5 = 1 should be positive. This 



3 



TABLE I: Eigenstates and energies of H\ + Hq F for two and three particles. Here Uj,Vj > 0, Uj + Vj = 1, « 2 = [1 — (A — 
J'/2)/n]/2, r 2 = [(A - J'/2) 2 + 2(J') 2 ] 1/2 , u§ = [1 - A/r 3 ]/2, and r 3 = [A 2 + (J') 2 ] 1/2 . A prime indicates new appearance of 
the same irreducible representation of the point group. States obtained by applying the spin lowering operator S~ or rotation 
of 7r/2 in the xy plane are not shown. 



notation eigenstate energy 
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implies, for two different t 2g orbitals in the atomic case, 
F - 5F 2 - 24F 4 > 0, or 

U - 3 J = U' - J > . (5) 

For a Slater determinant with both e g orbitals one ob- 
tains Fq — 8F2 — 9F4 > 0M> This condition is expected to 
be more restrictive than Eq. Q , since F 2 is usually more 



than one order of magnitude larger than F4. For exam- 
ple, using the above estimates for F 2 and F4 this gives 
F > 0.94 eV and V - J > 0.17 eV. We assume Eq. © 
to be valid in general. Otherwise for large |A| there is 
a charge-transfer instability of the ground state for the 
undoped system. The physics of Sr 2 Ru04 suggests that 
A is small and negatived 



III. EIGENSTATES OF THE LOCAL 
HAMILTONIAN 

The local part H\ + H l CF can be easily diagonalized. 
To describe the undoped system, we need the eigenstates 
with three electrons, and those with two and four elec- 
trons are needed when the effects of the hopping term 
Hh or doping are included. We denote the eigenstates 
by \inTSM), where i is the site index, n is the number 
of electrons, Y denotes the symmetry (irreducible rep- 
resentation of the point group D^h or symmetry of the 
basis function for the two-dimensional representation), S 
is the total spin and M its projection on the tetrago- 



nal axis z. If S = 0, M is suppressed. For simplicity, 
we drop the site index in this Section. The subscript g 
is also dropped in the irreducible representations. Some 
eigenstates and its energies are listed in Table I. The 
remaining ones for n = 3 and n = 2 are obtained by ap- 
plying the operator S~ or a rotation of tt/2 around z to 
those listed. The corresponding results for n = 4 can be 
obtained from those of n = 2 using electron-hole symme- 
try: replace creation by annihilation operators with the 
opposite spin, d\ aa — > di a ^- a , replace the vacuum by the 
state with all t 2g orbitals occupied, change the sign of A, 
and add U + 4U' — 2 J — 5U' to the resulting energies. 
Since we started with a local interaction Hamiltonian 



4 



Hj with full rotational symmetry, the symmetry group 
of H l j + H l CF is actually higher than D^. For example 



rotations around z of the orbitals yz and zx keeping the 
xy fixed leave H} + Hq F invariant. As a consequence 
the state |2&i0) is degenerate with 1 262O) , and the states 
|3aiiM) and \Za2hM') also become degenerate. This 
degeneracy is broken if the conditions U — U' + 2 J, J' = 
J are relaxed, but a degeneracy between \3ai^M) and 
\3bi^M') persists which is broken only if an exchange 
interaction J between the orbitals yz and zx different 
from the other two is introduced. 

For n — 3, the ground state is the spin quadruplet 
|3&i§M) if I A| < VT5J, while for | A| > Vl5J the ground 
state is also four-fold degenerate, but it is the spin and 
orbital E g doublet \2>x\M), \3y^M). These two possi- 
bilities lead to two different effective Hamiltonians in the 
undoped case, after integrating out the charge degrees of 
freedom. 



IV. EFFECTIVE HAMILTONIANS FOR 
UNDOPED PLANES 

In this Section we construct effective Hamiltonians H c s 
for the undoped case, using second-order degenerate per- 
turbation theory in Hh- Depending on the ground state 
of H — Hh, there are two possibilities for H e g. 



A. A < 15J 

In this case the ground-state manifold of H — Hh is the 
spin quadruplet |i36i|M) at each site i. The degeneracy 
is lifted by second-order contributions in which the in- 
termediate states have two nearest-neighbor sites with 2 
and 4 electrons, both with total spin S = 1 and both with 
the same symmetry B 2g , x or y. The different matrix el- 
ements are easily calculated using Eq. (0J and Table I. 
We omit the details. The resulting H c g is a Heisenberg 
model for the effective S = 3/2 spins: 



jja 
H cff 



K 



A(t 2 +t' 2 ) 
9(U + 2J) 



(6) 



The coupling constant K turns out to be indepen- 
dent of A. The ground state of this model is a 
two-sublattice antifcrromagnet with antiparallel nearest- 
neighboring spins. We call it AFI. The energy per site 
can be calculated accurately enough using spin waves, 
and is given by2& 



B. Large A 2 



H l CF is the 



For A 2 > 15J the ground state of H\ 
spin and orbital doublet \i3j^M), with 7 = x or y, which 
we denote briefly as \i^/u). The number of intermediate 
states is much larger than in the previous case, and H c g 
becomes very complicated. Since for A 2 > 15 J the struc- 
ture of the states involved in the derivation is already 
very similar to that for A — > ±00 (as can be checked by 
inspection of Table I), we restrict the calculation to this 
case. The result is 
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^ ^ P'ix Pi-\-S,y (JpSi.Si+s + A)} , (8) 

.5 

where now Si are spin-1/2 operators, are the orbital 
projectors 



(9) 



and 



Ja 
A 



2t 2 { 
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Jf 



2Jt 2 



U> 



J 2 ' 



(10) 



The first term of is a one-dimensional interaction, 
FM in the orbital degrees of freedom, and AF in spin. 
The term proportional to A is a spin independent AF 
orbital interaction, while the Jf term is FM in spin and 
AF in orbital variables. 

Clearly, there are two possible competing ground states 
of H^ s : i) a ferromagnetic orbital ordering (all sites \ixa), 
for example), with spin degrees of freedom determined 
by the critical one-dimensional AF Heisenberg model; ii) 
a spin FM and orbital Neel ordered phase (for example 
\ix 1) in one sublattice and \iy f) in the other). In the 
first case, for finite A, a small AF coupling between 
the chains appears [see Eqs. ((T3^l and (|T4f> ]. which yields 
long-range order at T — 0. We call this phase AFII. 
From Bethe ansatz results^ the energy of this phase for 
A — * ±00 is 



-Eafii = E ixl / 2 — Ja1ii2 



(11) 



The second phase will be denoted FM-AFO, and its 
ground-state energy is 



E 



FM-AFO 
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E 3xl/2 = -2(A+^) = - wr — J . (12) 



E A 



FI 



E 



3bi3/2 



-2KS(S + 0.158) = -4.974^, (7) 



THE PHASE DIAGRAM 



where E 3bl3 / 2 is the energy of the state |i36i|M) given 
in Table I. 



We now turn to the construction of a phase diagram 
for undoped planes, comparing the energies of the phases 
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U'/J 




3.2 3.4 3.6 3.8 4.0 4.2 4.4 

lAl/J 



FIG. 1: Phase diagram of undoped Ru02 planes for t — t' — 
J/4 (thin line) and t = t' = J (thick line). The phases are 
labeled as explained in the text. 



described in the previous Section, but now for arbitrary 
A. Since the correction terms for yl5J < | A| < +oo are 
small, we do not expect any new phases to appear in this 



interval, except perhaps near the borderline between two 
phases, as we will discuss at the end of this section. 

The energy of the AFI phase is still given by Eq. J7J, 
due to the independence of H£ s [Eq. (JHJ] on the crystal- 
field sppliting. The most important change occurs in 
the part of H^ s [Eq. that describes the AFII phase, 
in which an interchain coupling is generated. Thus, the 
effective Hamiltonian of the AFII phase, assuming that 
the chains are oriented along the x direction, is 



J^Si.Si+Sr + C} . (13) 



where J' A is the AF coupling along the chains, which 
coincides with J a for |A| — > oo, J± — 4(Cfj — Cu), 
C = Cf j + Cfi, and (Cjj_) is the negative correction 
of the energy due to virtual hoppings from site i to i + y 
and back when the spin of both sites are equal (opposite). 
The expressions for J' A , C aa i, as well as the energy of 
the FM-AFO phase are lengthy but straightforward to 
obtain, and we do not reproduce them here. Instead, Jj_ 
takes a simpler form, 



X (A 2 + J 2 ) [u 1 + 2(r 3 - J)] (C/' 2 + 4 [A 2 + J 2 + U'(2r 3 - J)]) ' 1 ' 

where r% = (A 2 + J 2 ) 1 / 2 . Eq. (J5J and the stability condition of the AFII phase against AFI for vanishing t and t' 
(A 2 > 15J) imply J±/J A < 6.80 x 1CT 4 . 



The energy of the AFII phase up to second order in 
Hh is given by 

E'afu = E 3xl/2 -J' A ln2 + C + E J± , (15) 

where Ej ± is the correction due to the interchain cou- 
pling. This correction can be calculated treating the J± 
term in a mean-field approximation, by a straightforward 
generalization of Schulz's results 38 for the case in which 
each chain has z nearest neighboring chains (he consid- 
ered z = 4 while we have z = 2). The energy gain due 
to the appearance of a spontaneous staggered expecta- 
tion value of the spin projection in the chain direction 
m = (— l) l (Sf) is (for any sign of J±) 

E J± = z\J ± \m 2 - 1(^)-V3 (z | J ± | m )4/3 . (16) 

Minimizing with respect to m one obtains the equilibrium 
value of the sublattice magnetization 

The resulting phase diagram is shown in Fig. ^ Due 
to electron-hole symmetry, the boundaries between the 



phases do not depend on the sign of A. The spin AF 
phase AFI, and the spin FM and orbital AF (FM-AFO) 
phase dominate the phase diagram. Comparing Eqs. i|ll|) 
and (|12f> . one obtains a critical value U' c ~ 8.52J for the 
boundary between the AFII phase and the FM-AFO in 
the limit |A| — * oo. For U' > U' c , the only stable phases 
are the spin AF ones. The main difference between these 
phases is that the staggered magnetization is very small 
in the AFII phase. In fact, from Eq. i|17|) we obtain 
m < 0.03 for the parameters of Fig. 2] 

Within our perturbation theory up to second order in t 
and t', the boundary between the FM-AFO phase and the 
AFII phase is independent of t and t' . The boundary of 
the AFI phase is also weakly dependent on hopping. For 
U'/J > 2, the stability region of the AFI phase is slightly 
enlarged by increasing the hopping parameters, as can be 
seen in Fig^ For smaller values of U'/J the energy of 
the FM-AFO phase decreases due to the proximity of a 
charge instability near which our perturbative treatment 
becomes invalid. The main effect of increasing t and t' 
is to enhance the energy difference between the stable 
and unstable phases in each region. These differences 
tend to be very small when the hopping parameters are 
small. For instance, the energies of the AFII and FM- 
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AFO phases completely coincide in the limit t = t' — > 
0. Thus, narrow-band systems are likely to show phase 
coexistence due to inhomogeneities. Also, we cannot rule 
out the appearance of more complex phases in a small 
region of parameters for which -Efm-Afo ~ -Eafii- One 
candidate is a phase in which orbitals display FM order 
in one direction (say x) and AF order in the y direction, 
while the spins are ordered antiferromagnetically in the 
x direction and ferromagnetically in the y direction. 



VI. THE DOPED SYSTEM 

Ru0 2 planes in RuSr 2 (Eu,Gd)Cu 2 8 are expected 
to have electron doping corresponding in our effective 
Hamiltonian H to a fraction below 20% of Ru sites with 
4 electrons. Depending on the ratio A/ J, there are three 
possibilities for the ground state of the local Hamiltonian 
H\ + H l CF for four electrons (see Table I) : 

1) A < — The ground state is the spin triplet 
\i±a 2 lM) (e.g., 1*4*211) = dl y] dl yi d\ yz] d\ zx] \Q)). 

2) < A < A c = (V41- 1) J/2 ~ 2.70 J — The ground 
state is the spin triplet and orbital doublet \i4xlM), 
\i4ylM). 

3) A > A c — The ground state is the spin singlet and 
orbital doublet |z4x0), |i4y0). 

Treating the hopping term in first-order degenerate 
perturbation theory, and combining with the results of 
Section IV, we can construct effective Hamiltonians H c g 
for the doped case. We begin by considering |A/J| < 
yT5 ~ 3.87, as suggested by the observed robust AF 
orders and the results of the previous Section. Then 
the ground state of H\ + H % CF for three electrons is the 
spin quadruplet with symmetry B\ g . 



A. -VTE.J < A < 

In this case, the problem of finding H c s reduces to cal- 
culating matrix elements of d[ di xya (all others vanish 
by symmetry) in the basis of |z3&i|Mi) and |£4a2lMo). 
For brevity we shall denote these states as |i|Mi) and 
|ilMo) . Using the Wigner-Eckart theorem, all matrix ele- 
ments can be calculated in terms of one of them (e.g., that 
for maximum projections, which is easily calculated), 
and Clebsch-Gordan coefficients ( JojMom\ J\Mi) for the 
combination of angular momenta Jo and j to give J\ . A 
similar approach was used before in problems of valence 
fluctuation with two magnetic configurations.- 9 Includ- 
ing the second order terms described before [Eq. ©], 
H e s becomes 

H% = -H E (l-Mo^MxKliM^lMO 

iS{M} 

X | (i + 8)^M[){{i + ^lMollilMoK^Mxl 



+Kj2(Si-Si-7)> (18) 

where {M} denotes the set M , Mi, M£, M[. Using the 
same method as above, it can be easily shown that this 
model is equivalent to a double-exchange model with in- 
finite antiferromagnetic coupling Jd e between localized 
and itinerant electrons: 

H °« = - f 'E (cUaH.c.) + j de ^a 

+ Kj2(Si-Sj-l)- (19) 

Here c\ a is the operator creating an itinerant electron of 

spin a at site i, and Sj = Yla3 c \ a <J aP c iP gi ves the spin 
of this electron. 

The physics of this model is expected to be quite sim- 
ilar to that of the model with FM exchange, as long 
as both Ru +5 and Ru +4 ground-state configurations are 
magnetic, which is this case. In fact, treating the spins 
classically as in Ref. [Fj, the sign of Jde is irrelevant for 
the electron dynamics, and only affects the effective mag- 
netic moment of Ru +4 . Thus, this results bring support 
to the model which successfully explained the magnetic 
properties of RuSr 2 (Eu,Gd)Cu20 8 «i^ 

B. < A < A c 

The 4d 4 configuration has orbital degeneracy in addi- 
tion to spin degeneracy. Proceeding as before, H e g can 
again be written as a double-exchange model, but now 
there are two types of carriers, each hopping only in one 
direction: 

ZtJ 

+J'^s,.S, + if^(S,.S J J). (20) 

' (ij) 

(2) 

While the magnetic properties of H^J should display 
some similarities to those of the previous H^§, anisotropic 
properties and the formation of stripes are more clearly 
expected here. No evidence of stripes in this system has 
been reported so far. A preferential direction was not 
observed in neutron experiments. 13 However, an equal 
amount of small domains with stripes oriented in the x 
and y direction cannot be completely ruled out by these 
experiments. Furthermore, although we explore here, for 
completeness, all the possibilities of our model, the avail- 
able experimental results seem to indicate that A < 0. 

C. A c < A < y/TEj 

In this case, the first order effective hopping vanishes. 
It is necessary to go to third order in Hh to effectively 
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exchange a 4d 3 spin quadruplet with a 4d 4 spin singlet 
and orbital doublet. Thus, the added electrons are es- 
sentially localized, and the observed magnetic properties 
would be difficult to explain within this picture. 

D. |A| > y/lEj 

Now, for any sign of A, the ground state of Hj + H % CF 
for n — 3 is the spin and orbital doublet \i3x^M), 
\iZy\M). We assume that the system is in the AFII 
phase to be consistent with neutron experiments, in 
spite of the different magnitude of the localized moment. 
Thus, the orbital degree of freedom is ferromagnetically 
frozen in the direction 7 = x or y. For negative A, 
H c g turns out to be equivalent to a double-exchange 
model with itinerant electrons coupled ferromagnetically 
to the localized spins 1/2, and one-dimensional hopping 
in the direction perpendicular to 7. Instead, for posi- 
tive A > V15J, the resulting H c g is equivalent to a t— J 
model with isotropic hopping and anisotropic exchange. 

VII. DISCUSSION 

We have studied the electronic structure of Ru0 2 
planes in RuSr 2 (Eu,Gd)Cu 2 0s using a strong-coupling 
approach to describe the 4d i 2ff orbitals of Ru and their 
interactions. For undoped planes (corresponding to for- 
mal valence +5 for Ru ions), we find three possible 
phases. Two of them are favored for large tetragonal 
crystal field (of any sign), and have orbital degrees of 
freedom which order at zero temperature (also at finite 
temperatures if hopping along the tetragonal axis were 
included). The spins order either ferromagnetically or 
in a particular AF order with very small staggered mag- 
netic moment compared to the experimentally observed 
one m ~ due to strong one-dimensional fluctu- 

ations. The dominant phase for small crystal-field split- 
ting consists of spins 3/2 which order antiferromagneti- 
cally, with nearest-neighboring spins pointing in opposite 
directions, as observed in neutron experiments^* One 
might wonder whether the effective measured staggered 
moment ~ 1.2 fig is closer to that of a localized spin 1/2 
rather than 3/2. However, for both mentioned AF phases 
there are several physical ingredients that reduce the mea- 
sured moment: i) spin fluctuations that reduce the sub- 
lattice magnetization, ii) effective Ru-Ru charge fluctua- 
tions (that are easily calculated within our perturbative 



approach), iii) Ru-0 charge fluctuations^ and iv) dop- 
ing, particularly if ferromagnetic polarons are formed^ 

When the system is doped with electrons, there are two 
main possibilities depending on the sign of the tetragonal 
crystal field parameter A. If it is negative (as it seems 
to be the case2* in Sr2Ru04), the additional carriers are 
described by a double-exchange model with infinite AF 
coupling with the localized 5 = 3/2 spins. This model is 
able to qualitatively explain the apparent contradiction 
between observed antiferromagnetic order, magnetic field 
dependence of the magnetization, and temperature de- 
pendence of the magnetic susceptibility^* It is also con- 
sistent with the observed magnetoresistanceiiS^*. Using 
previous results of the effective double exchange modelia, 
the experimental positive Curie constant O = 100 ± 3 
K suggests that the xy hopping t' ~ 0.25 eV. A more 
quantitative description of the magnetic properties re- 
quires an accurate calculation of the magnetic moment. 
It is also possible that the double-exchange model should 
be supplemented by inter-atomic Coulomb repulsions of 
a moderate ranged since the number of carriers in the 
system is low, particularly taking into account the low su- 
perconducting critical temperature. Also, previous stud- 
ies of the double-exchange model suggest that there is 
macroscopic phase separation at small doping 2140 which 
is inhibited by long-range Coulomb repulsion. 

If A is positive, the effective model for the doped 
case is similar, but the carriers have an orbital degree 
of freedom, as in manganites^i**^ which might lead to 
the observation of orbitons by Raman scattering** for 
enough doping. We also expect the formation of stripes in 
this case. However, there is no experimental evidence of 
stripes so far in the system, and fitting of optical proper- 
ties of another layered ruthenate Sr 2 RuC>4 suggests that 
A is small and negatived More detailed studies of the 
effects of doping would be possible if the superconduct- 
ing critical temperature T$ could be further enhanced, 
either by appropriate substitution of the rare earth or by 
applied pressure. 
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